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ABSTRACT 



The general second-order massive field equations for arbitrary positive integer spin in 
three spacetime dimensions, and their "self-dual" limit to first-order equations, are shown 
to be equivalent to gauge-invariant higher-derivative field equations. We recover most 
known equivalences for spins 1 and 2, and find some new ones. In particular, we find a 
non-unitary massive 3D gravity theory with a 5th order term obtained by contraction of 
the Ricci and Cotton tensors; this term is part of an j¥ = 2 super-invariant that includes 
the "extended Chern-Simons" term of 3D electrodynamics. We also find a new unitary 
6th order gauge theory for "self-dual" spin 3. 
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1 Introduction 

A peculiar feature of relativistic gauge field theory in three spacetime dimensions (3D) is 
that the particles associated with the gauge field may be massive even if the gauge 
symmetry is unbroken. In the well-known examples, such as "topologically massive 
electrodynamics" (TME) [IHI] or "topologically massive gravity" (TMG) [5], the non- 
zero mass arises because of some parity-breaking topological term in the action, but the 
recent example of "new massive gravity" (NMG) [6] demonstrates that the phenomenon is 
more general than had previously been thought; NMG is a generally covariant and parity 
preserving theory for interacting massive spin 2 particles with two polarization states of 
helicities ±2. An obvious question is whether there are other massive 3D gauge theories 
waiting to be discovered, and whether there is some systematic way of finding them. The 
main purpose of this paper is to show how gauge theories describing massive particles of 
integer spin may be found in a systematic way starting from standard second-order (or 
sometimes first-order) free field equations. Along the way we will recover a number of 
known gauge theory models, but we will also find a few new ones of interest. 

Let us first recall the Fierz-Pauli (FP) field equations [7] for a symmetric rank-s 
tensor field describing massive particles of integer spin s on a four- dimensional Minkowski 
spacetime (with 'mostly plus' signature): 

(□ - m 2 ) m ... Ms = , c^.-*-! = , ?T<W-p s - 2 = . (1.1) 

The differential subsidiary condition is needed for positivity of the energy (and hence 
unitarity) and the tracelessness condition is needed to avoid the propagation of lower 
spin modes; see e.g. [8]. The same equations may be used if the Minkowski spacetime is 
3-dimensional, and in this case they may be rewritten as 

\&{m)2?{-m% x "4> pn = , ?T<W-p s - 2 = , (1-2) 
where the operator 
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is an on-shell projection operator: 



^ 2 ( m )0 = &>(m)(j> if (j) satisfies (O]) . (1.4) 

Note that the equations (11.21) imply the differential subsidiary condition, and that the 
tensor ^(m)^(— m)0 is symmetric on its s indices as a consequence of the subsidiary 
conditions. The operator projects onto modes of helicity ±1 within the space of 
divergence-free vector fields satisfying the mass m wave equation (the sign depending on 
the sign of m) and hence onto modes of helicity ±s when acting on the rank-s symmetric 
tensor <fi. The equations (II. 2p therefore propagate one mode of helicity s and another of 
helicity — s, both with mass m. These are the two helicity states of a massive particle of 
spin s in three spacetime dimensions. Actually, each helicity state yields an irreducible 
induced representation of the 3D Poincare group but a parity preserving theory of spin 
s must propagate both helicities ±s with the same mass because parity flips the sign of 
the helicity. 

There is a generalization of the 3D FP equations (II .ip that breaks parity. Consider 
the following equations for independent positive masses m±: 

[^(m + )^(-m_)] w »<f> pn = , tT<W.. Ps _ 2 = . (1.5) 

These equations propagate one mode of helicity s with mass m + and another of helicity 
—s with mass m_; we shall refer to them as the "generalized (3D) Fierz-Pauli" equations. 
This generalization changes the dynamical equation from the mass m wave equation to 
the new equation 

(□ - m 2 ) w ... Ms = fie^ TU d T (t) u ^...^ , (1.6) 

where 

m 2 = m + m_ , fx = m_ — m + . (1-7) 

The subsidiary conditions are unchanged (as long as m 2 ^ 0) and these imply that the 
right hand side of (II. 6p is symmetric in its s free indices. We may assume without loss 
of generality that ft > 0, and take the limit that \x — > oo for fixed 

fi = m 2 /ft, (1.8) 

which implies that m 2 — > oo. Equivalently, we take m_ — > oo for fixed m + . In this 
limit the helicity — s mode decouples and we have a single mode of helicity s and mass /i 
described by the first-order equation 

^(/i) M /0 PAt2 ... Ms = , rT<W. p s - 2 = • (1-9) 

The differential subsidiary condition is again implied. Following the terminology intro- 
duced for the spin 1 case in [H] and subsequently used for the spin 2 case in [TU] , we shall 
call this the "self-dual" theory of a spin s particle of mass [i in 3D; the arbitrary spin 
case was studied in [TT] . 

We will show how these various sets of equations for a rank-s tensor may be replaced, 
in two different ways when s > 2, by equations for a rank-s gauge potential. The method 
was first used for s = 2 in [12] and we will shortly summarize the results obtained there. 
It also applies for s — 1, in which case the FP equations become the Proca equations; 
we shall use the notation B in place of <fr for the Proca vector field. Although this case 
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is too simple to fully illustrate the method, it is a good starting point. Specifically, let 
us start from the self-dual spin 1 model, for which the (unitary) action is 

Ssd[B] = J d 3 x {\e^B p d v B p + \^B 2 } . (1.10) 

One may easily verify that the field equation is 

^(/i)/^ = 0, (1.11) 

and that this implies the differential subsidiary condition d ■ B = 0, which has the 
following general solution in terms of a gauge potential A, defined up to a Maxwell gauge 
transformation: 

B» = \e^»F vp = F" , F pv ee 2d {il A v] . (1.12) 
When this is used in (11. lip , we find the new equation^] 

^(/i)/F„ = 0, (1.13) 

which is the equation of motion of TME: 

Stme[A] = J d 3 x [\F 2 + liie^A^A^ . (1.14) 

This establishes on-shell equivalence of the self-dual and TME theories; i.e. the equiv- 
alence of their field equations. Off-shell equivalence then follows because we are free to 
choose the sign (and have done so) such that the TME action is also unitary; in other 
words, the one propagated mode is physical, rather than a ghost. 

We should stress that the equivalence of the self-dual and TME theories is not a new 
result; the off-shell equivalence was proved in [13J by means of a "master action". It was 
also observed in [13] that the equations of TME are those of the self-dual model with 
A — > F(A), but the on-shell equivalence does not follow from this fact alone because the 
same replacement applied to TME yields the inequivalent equations of "extended topo- 
logically massive electrodynamics" (ETME) [T4], as we review in section |5J In contrast, 
the procedure of solving subsidiary conditions manifestly yields equivalent equations, al- 
though there is no guarantee that the corresponding actions are equivalent because it 
may be that the new action is not unitary. This point is well illustrated by the the 3D 
Proca theory for a spin 1 particle of mass m: the field equations are 

(□ - m 2 ) B p = , d-B = 0. (1.15) 

Solving the subsidiary condition as before, we arrive at the "extended Proca" (EP) 
equation (□ — m 2 ) F^ = 0, which is derivable from the following action containing the 
"extended Chern-Simons" (ECS) term introduced in [TJ] 

Sep[A] = f d 3 x [\e^ p F,d u F p - \m 2 e^ p A^d u A p } . (1.16) 

By construction, this model is on-shell equivalent to Proca, and this implies the prop- 
agation of two modes of mass m, one with helicity +1 and the other with helicity — 1. 

1 It is not generally permissible to substitute the solution of a field equation into the action, and no 
such substitution will be contemplated in this paper. 
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However, the on-shell equivalence to Proca does not extend to an off-shell equivalence 
because one of the two helicity modes propagated by the EP theory is a ghost, as we show 
later. We believe that this clarifies the very brief discussions of this model in [EUE)], 
where it was merely observed that it leads to "massive propagation of the field strength" . 

Moving on to spin 2, let us apply the analogous 'trick' to the self-dual spin 2 model 
with field equations 

^)/</v = o, ?T</v = o. (1.17) 

In other words, we solve the differential subsidiary condition d^4>^ u = that these equa- 
tions imply. The general solution involves a second-order differential operator because 
the symmetry of M „ means that it is divergence-free in both indices if it is divergence-free 
in one (as the subsidiary condition says it is). In fact, the general solution is 

for some symmetric tensor potential h^. The matrix operator £f is the "Einstein opera- 
tor" in the form introduced in [6]; the terminology comes from the fact that G>i n ' ) is the 
linearization of the Einstein tensor for the metric g^ v = r]^ + h^ v . The self-dual field 
equation of (11. 1TH now becomes 

^(/VGjfO = I) , ^ G (Hn) =0; (L19) 

but the second of these is implied by the first, which is the linearized limit of the field 
equation of "topologically massive gravity" (TMG): 

G^ v -\ — C^ v = 0, ■\f\g\C IJtV = e^ Tp D\ T G p \ v ) , (1.20) 
/i 

where \g\ — — det g for metric tensor g. The tensor G is the Einstein tensor and C is the 
Cotton tensor, the 3D analog of the 4D Weyl tensor. The TMG action is the sum of a 
Lorentz Chern-Simons (LCS) term, the variation of which is proportional to the Cotton 
tensor, and a 'wrong sign' Einstein-Hilbert (EH) term (i.e. with the sign opposite to the 
usual one; the sign itself is convention dependent). TMG is unitary despite the fact that 
the LCS term is 3rd order in derivatives because one may choose the overall sign of the 
action such that the one propagating mode is physical, rather than a ghost; this fixes the 
sign of the EH term. We should again stress that the equivalence of linearized TMG to 
the self-dual spin 2 theory is not a new result. It was proved by a 'master action' method 
in [HI], where the first-order self-dual spin 2 model was introducecll. 

We can repeat the above procedure for the FP spin 2 theory itself, and its 'generalized' 
version, and this was done in [12]. The resulting 4th order equations are the linearized 
equations of "new massive gravity" (NMG) and its extension to "generalized massive 
gravity" (GMG) ®. NMG has the action 

K = R^R^-Ir 2 , (1.21) 



2 Linearized TMG was also proved to be equivalent to a second-order self-dual model in [TBJ; this 
action is the m 2 — > limit of the action (13.11) to be considered below. The equivalence of all three spin 
2 models was made manifest in |17j . and a simplified version of this 'triple' equivalence was presented 
in [6j. For an earlier overview of spin-2 equivalences, with an emphasis on implications for interactions, 
see [32]. 



I NMG OC / d X V|# 
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where R pv is the Ricci tensor and R the Ricci scalar. GMG has, additionally, an LCS 
term with coefficient so one recovers NMG in the limit fi 2 — > oo for fixed m and 
TMG in the limit m 2 — > oo for fixed \i. For NMG, this proof of linearized equivalence 
with FP is closely analogous to the proof presented above that extended Proca is on- 
shell equivalent to Proca. That analogy, and the fact that we now have 4th order field 
equations, might lead one to expect one of the two spin 2 modes of NMG to be a 
ghost. Remarkably however, NMG is unitary. This was originally shown by an auxiliary 
field method that allows one to pass directly from the linearized NMG action to the 
unitary FP action [6|IT9] ; the absence of ghosts was subsequently confirmed by a canonical 
analysis |20j . The auxiliary field method does not seem to be as useful in the GMG case 
but the absence of ghosts in this case was demonstrated by canonical methods in [T2] . 

The method that we have now illustrated for spin 1 and spin 2 can be extended 
further in the spin 2 case [12J: we may solve both the differential subsidiary condition 
and the algebraic tracelessness condition simultaneously by writing 

</V = e { ; x d\ T ^ x \ v) ^h pu = Cf^ , (1.22) 

where the symmetric tensor potential h is now defined up to a linearized diffeomorphism 
and a linearized Weyl transformation. If this is used in the self-dual equations (I1.17p . 
one finds 4th order equations that are the linearization of the equations obtained from 
GMG by omitting the EH term. This model was called "New Topologically Massive 
Gravity" (NTMG) in [12] ; by construction, it is equivalent to the spin 2 self-dual model, 
as a canonical analysis confirms. The same model was studied independently in [2~T] . 

Solving both subsidiary conditions for the generalized FP theory yields 5th order 
equations that also have an obvious non-linear generalization; the 5th order term in 
the action is obtained from the contraction of the Ricci tensor with the Cotton tensor. 
Another purpose of this paper is to investigate this new spin 2 field theory. Let us recall 
here that the Cotton tensor C pv can be be written in terms of the 3D Schouten tensor 
S pu as follows 

Cfjn, = — .Ejj, p D T Sp U , S pu = R pv -^g pv R. (1.23) 

Thus, we find ourselves having to consider a contribution to the Lagrangian density of 
the fornjfl 

&elcs = VlglR^C^ EE S^S°D v S pa . (1.24) 

The Schouten tensor has an interpretation as a gauge potential for "conformal boosts" , 
so we see that the new 5th order term is a type of Chern-Simons term. We shall call it 
the "extended Lorentz Chern-Simons" (ELCS) term because it is analogous to the ECS 
term of electrodynamics. Indeed, we find that the ELCS term leads to ghosts in exactly 
the same way as does the ECS term. 

This result is nicely explained by the jV = 2 supersymmetric extension of 3D gravity 
models with an ELCS term. The JV = 2 supersymmetric extension of linearized GMG 
was presented in [12]; it unifies this 4th order 3D gravity theory with the 2nd order theory 
for massive spin 1. Here we show that the inclusion of the ELCS term in the gravity 
sector implies the inclusion of the ECS term in the electrodynamics sector, because both 



3 We understand that models including this term have been investigated in unpublished work of S. dc 
Haro and R. Jackiw. 
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are part of a single j¥ = 2 off-shell supersymmetric invariant. In addition, we obtain in 
this way a similar result for spin 3/2; one that is of course already implicit in the j¥ = 1 
supersymmetric extension: the inclusion of the ELCS term (ll.24p implies the addition 
to the Lagrangian density of a 4th order term of the form 

JS?=VW«^> (1-25) 

where if is the 'gamma-traceless' Cottino vector-spinor super-partner of the Cotton 
tensor [221123], and ^ M is its Majorana conjugate (we suppress spinor indices); upon 
linearization |12j . 

%n) = rd^ lin) + e^d u Mf^ , J>g in) = er'M, , (1.26) 

where ip^ is the anti-commuting Majorana Rarita-Schwinger (RS) vector-spinor potential. 
Obviously, the term (I1.25P must lead to spinorial ghosts because they are needed to 
complete the on-shell ghost particle JV = 1 supermultiplets implied by the existence of 
spin 2 ghosts. 

Although it is far from clear how one might introduce interactions consistent with 
higher-spin gauge invariances, it is of interest to consider how the procedure that we have 
been outlining for spins 1 and 2 generalizes to higher spins, and one purpose of this paper 
is to present some details of the application to spin 3. Obvious questions are whether 
there are analogs of the linearized TMG, NMG and NTMG actions for spin 2. A spin-3 
analog of TMG was proposed in [21] but this was based on the assumption (standard 
for 4D [25]) that the symmetric-tensor parameter of the spin-3 gauge transformation 
should be traceless. This allows the construction of a 2nd order analog of the linearized 
Einstein tensor, and hence a 3rd order TMG-like action, but this action was shown in [21] 
to propagate a spin-1 ghost in addition to the desired spin-3 mode. Our method leads 
naturally to a 3rd order analog of the Einstein tensor, and hence to a 4th order analog of 
TMG; although we present only the field equations of this topologically-massive spin-3 
theory (the action likely requires auxiliary fields) we confirm that it propagates a single 
spin-3 mode. There is also a natural spin-3 analog of linearized NMG but, as in the 
spin-1 EP case, it propagates one of the two spin-3 modes as a ghost. Finally, we derive 
the spin-3 analog of linearized NTMG, off-shell as well as on-shell. This is a 6th order 
field theory that is both unitary, as we confirm by a canonical analysis, and propagates 
a single spin-3 mode. 

The organization of this paper is straightforward: we proceed sequentially through 
spins s = 1, 2, 3. There is not much left to say about spin 1, but we complete the job in 
section |2J In section [3] we focus on the new 5th order spin 2 theory, and we show how 
jV = 2 supergravity combines the 3rd order ECS term of 3D electrodynamics with the 
linearization of the 5th order ELCS term for spin 2. We then adapt the procedure to spin 
3 in section HJ this involves the introduction of 3D analogs of the Einstein and Cotton 
tensors. We summarize in section [5] and include a brief discussion of higher spins and 
other possible generalizations. 
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2 Higher-derivative 3D electrodynamics 



We will now complete our discussion of spin 1 models by considering the generalized 
Proca theory with action 

S[B] = J d 3 x {±G 2 - Ifie^B^Bp - \m 2 B 2 } , G M = e^ p d v B p . (2.1) 

The field equation is 

[^(m^^i-m-)]/ B v = , m 2 = m+m_ , p, = m_ - m + > . (2.2) 

This propagates the helicity +1 mode with mass m + and the helicity —1 mode with mass 
m_. One may verify that the Hamiltonian density is positive definite, so these modes 
are propagated unitarily. Taking m_ — > oo for fixed m + (or jl — > oo for fixed fi = m 2 /jl) 
yields the self-dual model, while setting m + = m_, i.e. jl — 0, leads to the parity- 
preserving Proca model. We have already seen how these special cases are equivalent to, 
respectively, TME and the non- unitary EP theory of ( II. 16ft . We now consider the generic 
case. 

Solving the subsidiary condition d ■ B = as in (I1.12p . we arrive at the equation 

[^(m + )^(-m_)]^F u = 0, F» = e^ p d v A p . (2.3) 
This equation follows from the action 

S[A] = ~ J (fx^-e^F^Fp + jiF 2 + m 2 e^ A^A^ . (2.4) 

This is the generic "three-term" electrodynamics model considered in [H]. We have now 
established the on-shell equivalence of this model to ( 12.1 ft and this gives us the helicity 
content, which was not worked out in [T3]. However, there is no guarantee of off-shell 
equivalence. 

Following [H] , but choosing the gauge diAi = 0, we may decompose the vector 
potential as follows: 

A = , At = -^L= 84 , di = e ii d j . (2.5) 

The action ( 12. 41) now reduces to 

S[ip,<p} = J d 3 x ^-ipD(p + ^(^ 2 + 0D0) +m 2 0^| . (2.6) 

The TME limit is one in which jl — » oo for fixed \x = m 2 /fl, which implies that m 2 — » oo 
too. In this limit, ip is auxiliary and can be eliminated to yield a unitary action for 
that propagates a single mode of mass \i. Setting jl — gives us the EP model that is 
on-shell equivalent to Proca; this model should therefore propagate two modes of mass 
m, and this is clearly the case. However, it follows from the obvious diagonalization of 
the action when jl = that one mode is physical and the other a ghost. 
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When /i is neither zero or infinity, we may define new variables (U,V) by 

^0 = 11 + V, ip =v ^V, (2.7) 

to arrive at the action 

S[U, V] = jd z x{\ (UDU - VUV) + m 2 (U + V) V + \fi 2 V 2 } . (2.8) 

We see that there is a ghost, in agreement with [14J. The Hamiltonian can be diagonalized 
and the ghost decoupled by taking its mass to infinity, but this just leads back to TME. 
Thus, the addition of an ECS term to a TME action, which by itself propagates a 
massive mode of helicity ±1, leads to a modification of the mass of this TME mode and 
the propagation of an additional massive ghost mode of helicity =Fl. 

The focus of [H] was the "extended topologically massive electrodynamics" (ETME) 
model with action 

Setme[A] = \jd 3 x |f 2 - I^F M ^F P } . (2.9) 

The field equation is 

□F" = fxe^ p d v F p . (2.10) 
Compare this to the TME field equation 

nA» = fie^ p d v A p + d»(d-A) . (2.11) 

We see that the ETME equation is obtained from the TME equation by replacing the 
gauge potential by its field strength. One might think that this implies their equivalence 
since a zero potential implies a zero field strength, and a zero field strength implies a 
pure-gauge potential, which is effectively zero in the context of a gauge theory such as 
TME. However, they are not equivalent: by setting m 2 = in ( 12. 8p . one sees that ETME 
propagates a massless mode in addition to a massive one [H]. The overall sign is such 
that the massive mode is a ghost but we could change the sign to make this physical, in 
which case this 'reversed' ETME would differ from TME by the additional propagation 
of a massless ghost. What we wish to emphasize here is that the procedure of replacing 
a field by its field strength is not one that invariably leads to equivalent field equations. 
There is no contradiction with our procedure of solving subsidiary conditions for the 
simple reason that TME has no subsidiary conditions. 

3 Higher-derivative 3D gravity 

For spin 2, we start by considering the following action for a tensor field e pv of no definite 
symmetry: 

S[e] =-~J d 3 x je^ W + ^ Tp e^d T e pu + ^ (e^e pu - e 2 ) J , (3.1) 



8 



where e = rj^e^, and £f is the Einstein operator defined in (j!.18p . Note that £f is, by 
definition, symmetric in both index pairs; as a consequence, only the symmetric part of 
e M1 , enters into the first term in (13. ip . If we write 

where is symmetric, then one finds that the field equations are equivalent to 

[^(m+)^(-m_)]/</v = 0, rT</v = 0, ^ = . (3.3) 



The 'auxiliary' vector t appears irrelevant but is needed for the action, as shown in [10] 
in the context of the self-dual case obtained in the limit m_ — > oo for fixed m + . One may 
verify that the Hamiltonian density associated to the action is positive definite, which 
implies that the two spin 2 modes are propagated unitarily. 

As explained in the introduction, the differential subsidiary constraint d^cj)^ = 
implied by the equations (13. 3p (assuming m? ^ 0) has the general solution (11.181) in 
terms of a symmetric tensor potential h. Making this substitution in the generalized FP 
equations (13. 3p . we arrive at the equations 

[^(m+)^(-m_)]/G^ n) (/i) = 0, R^' m \h) = 0. (3.4) 

These are precisely the field equations of the linearized GMG model obtained by adding 
an LCS term to the parity-preserving NMG model. This guarantees on-shell equivalence 
of GMG to the generalized FP theory. It does not guarantee off-shell equivalence, i.e. 
the absence of ghosts, but this was shown in [12] following the canonical analysis used 
for NMG in EHI. 



As also explained in the introduction, one may solve both the differential and the 
algebraic tracelessness condition by equating <p to the linearized Cotton tensor for a 
symmetric tensor gauge potential h that is now defined up to the combination of a 
linearized diffeomorphism and a linearized Weyl transformation. In this way the first- 
order self-dual spin 2 equations become the linearized equations of the 4th order NTMG 
theory [12] . Here we shall use the same procedure to convert the second-order generalized 
FP equations for spin 2 into equivalent 5th order equations. Specifically, we use (11.221) 
in the 2nd order equations (13.31) to arrive at the equivalent 5th order equation 

[^(m+)^(-m_)]/C^ = 0, (3.5) 

which is invariant under the gauge transformations 

hpv -> + d(jj,v v ) + T)nv u (3.6) 

for arbitrary vector field v and scalar field u. The subsidiary conditions have now been 
replaced by the gauge-invariant identities 

&>C<M = , rTCff> = . (3.7) 
The gauge-invariant action is 

S[h] =nj d 3 x {- Vc^ - ^e^d T C^ + ±hn3C<Sp} . (3.8) 
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' J ^r,<^- h 00 = —N, (3.11) 



This is the quadratic approximation to the 5th order action 

S = fi J d 3 x {-J?lcs + ^VW\K - ^VWlR^cA , (3.9) 

which includes (for finite m 2 ) an ELCS term. Note that the linearized Weyl invariance 
of the linear theory does not generalize to the full theory. Instead one finds a conformal 
covariance property, explained for NMG in [6]: a scalar density is conformally covariant 
if the Weyl variation of its integral is the integral of a function times the same scalar 
density. This is a property of the scalar density y\g\ K of NMG and it is also a property 
of the ELCS term. As both these scalar densities have no linear term, their conformal 
covariance implies a linearized Weyl invariance of the quadratic approximation. 

We now have a new 5th order theory for spin 2 but there is no guarantee that neither 
spin 2 mode is a ghost. We shall investigate this issue in the context of the linearized 
theory with quadratic action (13. 8p . The six components of the metric perturbation 
are determined up to the gauge transformation (13. 6p by two functions. We may choose 
the gauge such that 

fcy = 0, $/km = 0, (3.10) 

which means that 

for two functions (£, iV). This decomposition can be obtained from that given in [12] 
(which is a variation of that used in |20j) by imposing the additional Weyl gauge condition 
S^hij = 0, and the results for the Weyl-invariant Cotton tensor components are therefore 
identical to those given in [12]. Here we go directly to the quadratic action, which 
becomes 

S[t,N) = lfJ d3x {-Z N + l(^ + \ N2 )+^ Nn t} ■ ( 3 - 12 ) 

In the m 2 — > oo limit for fixed \x we recover the result of p2] for the NTMG model; the 
variable N is auxiliary in this limit and its elimination yields an action that propagates 
unitarily a single mode of mass //. 

If we first divide by the overall factor of fi in (I3.8P then we may take the fi — > oo limit 
for fixed m 2 to get a 5th order model with the following Lagrangian: 

l = -\hrc^ + ^rncjjo . (3.13) 

By construction, this is on-shell equivalent to the standard 3D FP model for massive spin 
2. However, applying to (I3.12p the same procedure of dividing by fi and then taking the 
H — > 00 limit we arrive at the 'reduced' Lagrangian 

L = ^- 2 N(n-m 2 )Z, (3.14) 

which shows that one mode is a ghost. The 5th order model with Lagrangian (I3.13P is 
therefore a much closer spin 2 analog of the EP model than is linearized NMG, so it could 
be viewed as the linearized version of an "extended Fierz-Pauli" theory. Alternatively, 
the non-linear theory with action 

S = fd 3 x [% LCS + ^y/^BTcA , (3.15) 
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could be viewed as "extended NMG" (ENMG) theory. 

Returning to ( 13 . 1 2 j) and assuming that /x is neither zero nor infinity, we may define 
new variables (U, V) by 

(=±{U + V), N = V2fxV, (3.16) 

to arrive at the equivalent action 

S[U, V] = Jd z x{\ (UDU - VDV) + m 2 (U + V) V + \fi 2 V 2 } . (3.17) 

This is identical to f l2.8|) . which means that the ELCS term leads to ghosts in 3D gravity 
in precisely the same way that the ECS term leads to ghosts in 3D electrodynamics. We 
next show how this fact is made manifest in jY = 2 3D supergravity. 



3.1 Higher-derivative jY = 2 supergravity 

We now aim to embed the action (13.81) into one of linearized jV = 2 3D supergravity. 
Because of the linearized Weyl invariance, we need only consider the Weyl supermultiplet, 
which consists of the fields {/i^,^,^}, where the (anti-commuting) RS fields ip® (a = 
1, 2) are Majorana spinors. These fields are subject to the following gauge invariances 
that reduce the number of independent off-shell degrees of freedom to 4 bosons and 4 
fermions: 

<5V = 9 (A)+V W > = 9 a*C° + , 5A )1 = d fl a. (3.18) 

The vector field is the parameter of linearized diffeomorphisms, the scalar fields u 
and a are the Weyl and abelian gauge parameters, respectively, and the Majorana spinor 
fields 7] a and ( a are the parameters of linearized S-supersymmetry and Q-supersymmetry, 
respectively. In the linearized theory, the Q-supersymmetry gauge invariance is indepen- 
dent from the rigid supersymmetry transformations, which are 



Sift = -\d p K^ a t a -\A»t ah e\ (3.19) 



This result agrees with [12] when account is taken of the linearized gauge invariances. The 
vector-spinors (ft® are the (dependent) S-gauge fields; the RS field-strengths are defined 
as in fll.26p . The transformations ( 13 . 1 9[) imply the following JV = 2 supersymmetry 
transformations of the conformal field strengths: 

^(Z) = l,e a C^ n) + le^e ab lalu e b d p F% (3.20) 
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Here, is the doublet of "Cottino" vector-spinors; as defined for jY = 1 in ( ll.26p . 

Using these transformation laws, one may verify invariance of the action 

S-^ 2 = J d 3 x [\hrC<M + |^ n a ) - \A»F») . (3.21) 

To check the invariance one needs the various Bianchi identities satisfied by the super- 
conformal field-strengths, the tracelessness of the Cotton tensor and the 'gamma trace- 
lessness' of the Cottino tensors. In agreement with [12], we see that the CS term of 3D 
electrodynamics and the LCS term of 3D gravity appear in the same JY = 2 invariant. 
However, the formula ( 13.211) can now be used to construct a whole family of JY = 2 
supersymmetric actions. 

For each set of fields that transform as in (13.191) modulo linearized gauge transforma- 
tions of the form (I3.18p . the action formula (I3.2ip yields an jY = 2 super- invariant. For 
example, one such set of fields is 

= e { ^d\ p h aW) , ^ W = e, vp d v r p , 4 1} = ^ UP duA p . (3.22) 

To be more precise, these derived objects transform as in ( 13.191) provided one ignores 
'pure- gauge' terms of the form ( 13. 18ft (which is justified because they automatically drop 
out of a gauge-invariant action). To verify this, we first note that ^ transforms into 
in the required way; the candidate for h$ found this way is not manifestly symmetric 
but this can be remedied by adding field-dependent Q- and S-gauge transformations. 
Similarly one finds that h$ does not transform into ip p ^' but instead that 

$h$ = \e a i, {i T d T r v + erd p r a ) + (3.23) 

However, the first term can be rewritten as the second term up to a Q- and S-gauge 
transformation by making use of the identity 

ef = ' 2^ [P 7 <T] • (3.24) 

Using (I3.22p in the action formula (13.211) yields the supersymmetric invariant 

sfr 2 = J d 3 x [\h^cp + m (1 \Z) - 14 1} ^} ■ ( 3 - 25 ) 

This is the JY = 2 supersymmetrization of the curvature-squared K invariant of NMG, 
and is again in agreement^ with the result of [12J . 

At this stage we have merely confirmed the results of [12] for the JY = 2 supersym- 
metrization of the linearized LCS and curvature-squared terms. However, we may now 
go further by iterating (13.221) to get a new set of fields, again up to pure-gauge terms: 

h$ = -2J2jf> , % (2) = \V« ^ , A$ = e^d v F p . (3.26) 

The transformations of these fields are again of the required form, so we may use them 
in the action formula (I3.2ip to get a new JY = 2 invariant that includes both the ECS 
and the ELCS terms: 

S&= 2 = j d 3 x {-RfaCjM + W ll)a ^iZ) - ¥" vp FAF P ) ■ (3.27) 

4 To compare, one should again use that the expression ipt-^ = e n" p 9uipp occurring in the action 
(|3.25|> differs from the expression j T drtp^ occurring in (7.9) of [12] by a Q- and S-gauge transformation. 
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This result implies that the ECS and ELCS terms lead to ghosts in precisely the same 
way since any ghosts must fill out on-shell jY = 2 particle multiplets. 

It should now be clear that we can construct, by further iteration, any number of 
still higher-derivative JV = 2 invariants. For example, the following set of fields has the 
required transformations 

h® = ~2C^ , ^ = \Yd T % ^ , Af = OF, . (3.28) 

Substitution into the action formula yields an JV = 2 invariant containing the FDF 
invariant of electrodynamics and the square of the Cotton tensor. 



4 Spin 3 

Let us start from the first-order self-dual equations for spin 3: 

^V)/</W = , ^ = tT</v p = . (4.1) 

These imply the differential subsidiary condition d^(f)p Vp = 0, which has the following 
general solution in terms of a 3rd rank symmetric tensor potential h: 

^y.vp ^[ivp h a p~y = G fii/p , (4.2) 

where 

% yp ^ = —e^e^e^drd^dr, (4.3) 

is the spin 3 analog of the Einstein operator. The tensor G pup is a spin-3 analog of the 
linearized Einstein tensor; it is invariant under the spin 3 gauge transformation 

h^up -> h^p + d^Zvp) (4.4) 

for arbitrary symmetric tensor field ^. The self-dual spin 3 equations are therefore equiv- 
alent to the 4th order equations 

<n«)/GW = 0, Gp = V ^G, up = 0. (4.5) 

The first of these equations implies the second, but only because of the lack of symmetry 
in its three indices; an alternative equivalent set of equations is 

&(n)(p a G U p )a = 0, Gp = 0. (4.6) 

By construction, these equations propagate a single spin 3 mode (and we verify this 
below). They are the natural generalization to spin 3 of the linearized TMG equations 
for spin 2. We expect the action to involve auxiliary fields, which we will not attempt to 
find here. However, it will obviously be possible to choose the overall sign of the action 
so as to ensure unitarity, and hence a 4th order unitary theory describing a single spin-3 
mode. 

Before proceeding to other cases, we should mention that a 3rd order action for a 
symmetric tensor gauge field was proposed in |24j as a spin-3 analog of linearized TMG. 
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Since the symmetric tensor parameter of the spin-3 gauge transformation was restricted 
to be traceless, it was possible to find a 2nd order rank-3 invariant tensor analogous to 
the Einstein tensor, and the equation proposed in [24J is formally the same as the first of 
the equations (14. 6p but with our 3rd order rank-3 tensor G replaced by their 2nd order 
rank-3 tensor. The resulting equation is derivable from an action, but this (3rd order) 
action propagates a spin 1 ghost in addition to the desired massive spin 3 mode [24J. As 
stated above, we believe that the 4th order equations (14. 6 p are the natural generalization 
to spin 3 of the linearized TMG equations. There is a further generalization to spin s 
that we present in the following section. 

We now turn to the standard FP equations for spin 3 

(□-m 2 )</v P = 0, d^p Up = 0, r/"Vw = 0. (4.7) 

Solving the differential subsidiary constraint, as in (14.2j) . we arrive at the 5th order 
equations 

(n-m 2 )G„ vp = 0, Gp = 0. (4.8) 
One easily shows that these equations are equivalent to 

2 

C*iJ,vp — 3^ = , (4-9) 

lib 

where 

1 3 

C^p = -HG^p - -Q(^G P ) , Qp U = j^D - dpd u , (4.10) 

Z o 

which is identically traceless: r\ vp Cp Vp = 0. A crucial property of both G and C is that 
they take the form 

Gp V p ^pvp ^h a ij^ , Gp Up ^pup ^h a ij^ , (4.11) 

for self- adjoint matrix operators £f and ^ . This means that the equation (14. 9 p follows 
from variation of the action 



S = - I </ : \r //""'' 





C£ afil Off a/37 

^>pUp n <9pvp 



h afh . (4.12) 



This is the spin 3 analog of the linearized NMG action for spin 2 and the extended Proca 
action for spin 1. The construction does not guarantee the absence of ghosts, so we must 
address this question by other methods. We shall see that one of the two spin 3 modes is 
a ghost. For this reason, we will not bother to investigate the gauge-invariant formulation 
of the generalized FP equations for spin 3. 

The rank-3 tensor C defined in f)4.10p is the spin 3 analog of the linearized Cotton 
tensor; it is invariant under the spin 3 analog of the linearized Weyl transformation 

hpyp -)> hp Up + i]^ u u p) , (4.13) 

in addition to its invariance under the transformation f)4.4p . and it satisfies the identities 

d^C, U p = 0, V ^Cp Up = 0. (4.14) 
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An alternative expression is 

1 3 

C/_ wp = — £{^ ra £j J ^ d\ T d a \S p } a p , S^yp = G^yp — —rj(p U G p -) . (4-15) 

The rank-3 tensor S is the spin 3 analog of the Schouten tensor. Note the explicit 
symmetrization that is needed in this expression for C, in contrast to the analogous 
formula for spin 2 in ( ll.23p . The 5th order action 

iconf / j3 



Sr = Jd^Cp Vp (4.16) 

is the "conformal spin 3" action presented in |26j, which has no propagating modes. 

As in the spin 2 case, we may also solve both the differential and the algebraic sub- 
sidiary conditions by writing 

<P^ P = Cp Vp . (4.17) 
We thus find that the standard spin-3 equations are equivalent to the 7th order equation 

(□ - m 2 ) Cp Vp = . (4.18) 

This is derivable from the action 

S = ^Jd 3 x W vp (□ - m 2 ) Cp Vp . (4.19) 

We shall show below that this propagates one of the two spin-3 modes as a ghost, so this 
spin 3 model is analogous to the EP model for spin 1. 

The problem of ghosts should not arise in the self-dual limit, in which only one mode 
is propagated. In this case, we use (I4.17P in (14. ip to arrive at the equivalent equation 

^(/i)AW = 0. (4.20) 

This is a 6th order differential equation for the gauge potential h, invariant under the 
gauge transformation 

Shfu, p = d { pi vp) + ri (p „ujp) , (4.21) 

for arbitrary symmetric tensor parameter £ and vector parameter u. It should be appre- 
ciated that this gauge potential is not the same as the one in (14. 2P ; their dimensions are 
different. The equation (I4.2(jp is derivable from the 6th order gauge- invariant action 

S = ±Jd 3 x h^&^^C^p . (4.22) 

We shall verify below that this action propagates, unitarily, a single (spin 3) mode. This 
is the spin 3 analog of NTMG. 

4.1 Canonical analysis 

We may choose the gauge 

d l h lup = . (4.23) 
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In this gauge we may write 



1 



i 000 



1 



h, 



OOi 



-v 2 y 



(-v 2 ) 

where we recall that, by definition 



2 



:-v 2 p 



jdidjdk^ 



di = e l3 d j . 



Note the following useful identities 

didj = SijV 2 - didj , didi = . 
A computation of the components of the tensor G^ vp yields the result 



* \ 2 



^000 

Goij 



(-v 2 ) 

6 



G 



00 i 



6 



G 



6(-V 2 )^ 
1 



ijk 



di<p 2 + di<f> 3 

didjd k (po + 3(%i9j(9 fc )0i + Sd^djd^fa + didjd k dffa 



6(-V 2 )2 

From this result we find that the components of G p are 

f-V 2 ^ 



G 
G, 



6 



(□03 - 0! 



+ ft D03 - fa 



6(-V 2 )^ 

Using these results, one may show that the equations (14.61) are equivalent to 
0o = —mfa = m 2 fa = -m 3 fa , ( n - m 2 ) fa = , 

which confirms the propagation of a single mode of mass m. 

Moving on to the components of C^p, we similarly find the result 



Cooo 
Coij 

Gijk 
where 



16 
1 



-fa 



00/ 



48 



diipo + 3difa 



16(-V 2 )^ 
1 

48 (-V 2 )^ 



didjDipx + -dQdjrfo + didj fa 



(4.24) 



(4.25) 



(4.26) 



(4.27) 



(4.28) 



(4.29) 



(4.30) 



didjdkUfa + 9d (i djd k pfa + Sd^djdk^o + 3didjd k d?fa 



^O = 0O + 3n0 2 , fa = fa + -D(f) 3 . 



(4.31) 
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We see that C avp depends on only two linearly independent combinations of the four 
variables (0 O , <t>i, 02, 03)- 

To analyse the Lagrangian (I4.12p we need the following results: 

1 , 



bT p G 



VP 



fivp 



h<P2 - g<P0<?3 



1 4 

-^102 - g^O03 + ^02^03 , 



(4.32) 



We now find that 



G uvp Q C 'uvp 



Am 2 



^o-02 "01 _ 



4m 2 



-0 2 (D-m 2 )0 3 . (4.33) 



The first term is zero on using the (ipo, equations. This leaves the second term, 
which propagates two modes of mass m but with one of the two being a ghost, as stated 
earlier. We see that the on-shell equivalence to the spin-3 FP theory does not extend to 
an off-shell equivalence. 

We now turn to the higher- derivative spin 3 theory with action (I4.19p . The Lagrangian 
depends only on the variables (ipo^i)'- 



L = -— ip Q (□ - m 2 ) • 



(4.34) 



We thus confirm that two modes of mass m are propagated, but one is a ghost, as stated 
earlier. 



Finally, we consider the action (14. 22ft . We may write its Lagrangian as 



1 



G uvp J aup 



L = W vp 

To analyse this Lagrangian, we compute 



J pup — £(ju <9|tCa|^p) 



h^Jpvp = tooo + 3^i n^ooo = ^Vo + 7^°^ • 

4o lb 



(4.35) 



(4.36) 



The first equality is obtained by using the identities satisfied by C avp . Using this result 
we have 



Cpi/p J uvp 



1 



(4.37) 



|g (ip ~ 3/i^i) 2 + jg-Vi (□ - /* 2 ) V>i • 

Eliminating the auxiliary field ipo removes the first term. What is left is a unitary 
Lagrangian for ipx that propagates a single mode of mass /i. We know that this has spin 
3 so (I4.22p is a unitary 6th order gauge-invariant action for a single spin 3 mode. 



5 Discussion 



Relativistic field equations that are higher than second order in derivatives are typically 
consigned to the class of "higher-order" equations exhibiting unphysical behaviour, asso- 
ciated classically with negative kinetic energies and quantum mechanically with negative 
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norm states, i.e. ghosts. However, it has been clear for some time, from the example of 
the third-order "topologically massive gravity" (TMG) [5J, that three dimensional (3D) 
spacetime allows exceptions, and the recently discovered "new massive gravity" (NMG), 
and its generalization to "general massive gravity" (GMG) which includes both TMG 
and NMG as special cases, has shown that there are even fourth-order equations that 
are free of the bad features usually associated with "higher derivatives" [6] . These grav- 
itational theories are unusual in that they are generally covariant but describe massive 
spin 2 gravitons. 

It was shown in [12] that the linearized versions of these 'massive gravity' theories 
can be found systematically by solving the differential subsidiary condition in standard 
second (or first) order 'non-gauge' massive spin-2 field theories of Fierz-Pauli (FP) type, 
and a "new topologically massive gravity" (NTMG) was also found this way. Here we have 
developed further this method, starting with a systematic presentation of the 'generalized' 
FP model for massive spin s and its "self-dual" limit. We used the simpler application 
to spin 1 to illustrate the method, and this led us to a discussion of models related to the 
"extended topologically massive electrodynamics" (ETME) of [H]. Because the method 
only establishes equivalence of linear field equations, there is no guarantee of off-shell 
equivalence, which amounts to the same as unitarity. This has to be checked by some 
other method, and a convenient one is the canonical decomposition method used in [T4] 
to show that ETME has ghosts. This method was also used in [2D] to verify the unitarity 
of NMG and in [12J to prove unitarity of (super)NMG; we have made extensive further 
use of this method here. 

Application of these methods for spin 2 led us to a novel 5th order spin 2 theory 
with a natural extension to an interacting generally covariant massive gravity theory. 
The 5th order interaction term is found from the contraction of the Ricci and Cotton 
tensors of the 3-metric. This term can be viewed as a type of Chern-Simons (CS) term 
and as a gravitational analog of the "extended CS" (ECS) term of ETME. In fact, this 
is more than an analogy: we have shown that both ECS and ELCS terms are part of the 
same JV = 2 supersymmetric invariant, which also includes the 'square' of the linearized 
Cottino tensor-spinor. As any ghost modes must form particle supermultiplets in this 
context, both ECS and ELCS terms lead to ghosts in precisely the same way, as we have 
explicitly verified. Although it might be considered disappointing that the new model for 
spin 2 fails to be unitary, we prefer to see this result as an indication of the uniqueness 
of GMG; there is a limit to higher derivatives if one insists on unitarity and it seems as 
though this limit is saturated by GMG. 

In the spin 3 case, we found a 4th order generalization of the spin-2 equations of 
linearized TMG. Obviously, these differ from the 3rd order equations, with a weaker 
gauge invariance, that were proposed in [23] as a generalization of TMG; in particular, 
our 4th order equations actually do propagate only a single mode of spin 3. We expect 
the covariant action to require auxiliary fields but if an action exists it will be possible 
to ensure unitarity by a choice of sign. In contrast, we found that the natural spin-3 
generalization of linearized NMG, which is of 5th order, has ghosts. However, we were 
able to find a unitary gauge-invariant 6th order action for a single spin 3 mode that is 
equivalent to the first-order self-dual theory for spin 3. This is a spin 3 analog of the 
NTMG model for spin 2 found in [T2J [2T] ; as in that case one proceeds by finding the 
general simultaneous solution of both differential and algebraic subsidiary conditions. 



18 



For spin 2 this is done by setting the FP field equal to a linearized Cotton tensor for a 
symmetric tensor potential, and the spin 3 case involves a rank-3 tensor generalization of 
this tensor. It seems likely that there is a spin-s generalization of this construction: there 
is a spin-s generalization of the linearized Cotton tensor that involves 2s — 1 derivatives 
[26] . so we expect the self-dual spin-s theory to be equivalent to a conformal-type gauge 
theory of order 2s. It is then natural to conjecture that 2s is the maximum number of 
derivatives consistent with unitarity for spin s. 

For any free field theory that propagates a single mode, one may always choose the 
overall sign of the action (assuming that there is one) such that this one mode is physical 
rather than a ghost. For this reason, unitarity is not really an issue in applications of 
our method that start from a self-dual model. It is only when there is more than one 
propagating mode that there is a potential problem: if there are modes with opposite sign 
kinetic energy then unitarity is lost for any choice of overall sign. The simplest example 
is the "extended Proca" (EP) gauge theory which is equivalent to Proca on-shell but not 
off-shell because one of the two spin 1 modes it propagates is a ghost. One plausible 
explanation for this feature is that the EP vector potential has opposite parity to the 
vector field of the Proca theory. In fact, one may see by inspection that in all cases 
considered here, for arbitrary s, ghosts arise in a model propagating both helicities if and 
only if the parity of the FP non-gauge field differs from that of the gauge potential. This 
leads us to conjecture that this is a general feature. In particular, we conjecture that the 
gauge theory will be unitary whenever the parity of the gauge field is the same as the 
original FP field. If true, this conjecture would nicely explain the unitarity of GMG, and 
it could simplify the construction of unitary gauge theories for higher spin. 

It is straightforward to take the first step towards a higher-spin generalization of our 
results. The general solution of the spin-s differential subsidiary condition is 

The rank-s tensor G is invariant under the spin-s gauge transformation 

(5.2) 

for arbitrary rank-(s — 1) symmetric tensor £, and it satisfies the Bianchi-type identity 

PG vltl ... llm _ 1 = 0, (5.3) 

which replaces the differential subsidiary condition. For example, the spin-s self-dual 
equations for a single mode of spin s are equivalent to the equations 

^ (w/xi P Gpfj,2---fj, 3 = > Gpi— ti a -2 — V P G upfll ... fls _ 2 = (5-4) 
This equations generalize to spin s the equations of linearized TMG. 

The linearized NMG equations can be similarly generalized. For example, the FP 
equations for spin 4 are equivalent to the 6th order equations 

(□ - m 2 ) Gyvpc = , vTG^ = • (5.5) 

This case is of interest in light of our unitarity conjecture because the parity of h is the 
same as the parity of <f> for all even integer spins, and the spin 4 case is the first instance 
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beyond spin 2. Thus, we expect there to exist a unitary 6th order parity-preserving 
action for spin 4. More generally, we expect there to exist for any even positive integer s 
a unitary parity-preserving action of order s + 2 that propagates both helicity modes of 
spin s. The problem is that it is far from clear how to find an action for these equations 
when s > 2; we would guess that it exists but involves auxiliary fields. 

Other possible generalizations are to theories of all higher spins, integer or half integer, 
as considered for 3D in [26H28] . and to fractional spin along the lines of [291 130] ■ It is 
also possible that there is some generalization to higher dimension but 3D is special in 
that the index structure of the independent field is preserved in the process of solving 
subsidiary conditions. The 4D Proca case illustrates the point: solving the differential 
subsidiary condition leads to the equation (□ — m 2 ) = 0, where now = e^^dyApu 
for antisymmetric tensor potential A. There is no action for A alone that yields this 
equation because there is no longer a match between the index structure of the potential 
and its equation of motion. 
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